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C3 ■ Abstract 




We study the structure of the constrained minimizers of the Gates-Lebowitz-Penrose free- 
energy functional ^ 7 glp(w), non-local functional of a density field m(x), x £ 7~l, a d-dimensional 
torus of side length L. At low temperatures, .Fglp is not convex, and has two distinct global 
minimizers, corresponding to two equilibrium states. Here we constrain the average density 
L~ d J T m(x)dx to be a fixed value n between the densities in the two equilibrium states, but 
close to the low density equilibrium value. In this case, a "droplet" of the high density phase 
may or may not form in a background of the low density phase, depending on the values n and 
L. We determine the critical density for droplet formation, and the nature of the droplet, as a 
function of n and L. The relation between the free energy and the large deviations functional 
for a particle model with long-range Kac potentials, proven in some cases, and expected to be 
true in general, then provides information on the structure of typical microscopic configurations 
, of the Gibbs measure when the range of the Kac potential is large enough. 
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1.1 The free energy functional 



Let Tl denote the d-dimensional torus with edge length L. Let M. denote the set of measurable 
functions m from Tl to [—1,1]. Here, m € A4 is an order parameter field representing the local 
average magnetization in an Ising system on a lattice in T^, viewed on a mesoscopic scale in which 
the microscopic lattice structure is invisible. The underlying microscopic model has a well-known 
lattice gas interpretation, in which p := (1 + m)/2, which takes values in [0,1], is viewed as a 
particle density that is bounded above due to the microscopic constraint that at most one particle 
may occupy any lattice site. 
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The object of this paper is to investigate the minimizers, and near-minimizers, of certain meso- 
scopic free energy functionals on Ai under a constraint on the value of the average magnetization 
(or density) L~ d J^- m(x)dx. These free energy functionals (defined just below) arose in work of 
Lebowitz and Penrose [25] and Gates and Penrose |18| . and have been proven to be large deviations 
functionals for the underlying Ising systems, at least in certain cases [TJ[30]. Thus "near minimal 
free energy profiles" m correspond to "typical" coarse grained magnetization configurations for the 
Gibbs measure of the underlying Ising system. 

We investigate the nature of the (near-) minimal free energy profiles m for sub-critical temper- 
atures, so that there are two equilibrium values m± of the magnetization, and for values of n just 
slightly above m_. Depending on the amount of the excess magnetization n — m_, the surface 
tension between the two phases, the bulk compressibility, and on L, the excess magnetization will 
either be uniformly dispersed over the entire volume, or will partly aggregate into a "droplet" of 
the m+ phase in a sea of the m_ phase. This phenomena was studied rigorously on the microscopic 
scale for the two dimensional Ising model with nearest neighbor interactions by Biskup, Chayes 
and Kotecky [6]. 

Here we consider this same droplet formation problem, but in any dimension d, and for Ising 
systems with long range Kac potentials instead of nearest neighbor interactions. While [6] presents 
direct analysis of the microscopic system, our starting point is instead to analyze a variational 
problem for a mesoscopic free energy functional where a number of tools from the calculus of 
variations can be used. Moreover, the connection between the microscopic Gibbs measure and the 
mesoscopic free energy functional has already been thoroughly investigated [HEQ], and we can rely 
on this to make contact with the microscopic scale. Finally, the mesoscopic variational problem is 
of interest in its own right. 

The Gates-Lebowitz-Penrose (GLP) free energy functional .Fglp that we study has the form 

^GLp(m) := / f(m(x))dx - - [ [ m(x)J(x - y)m(y)dxdy , (1.1) 

where / is some convex function, and J is a non-negative radial function with 

/ J{x)dx = 1 . (1.2) 

We further suppose that J has finite range range, which we then use to set our length scale so that 
J has unit range. More specifically, we require J(x) = for |x| > 1, and J(x) > a, some a > 0, for 
\x\ < 1/2. 

Such functionals were introduced by Gates and Penrose [18], building on previous work by 
Lebowitz and Penrose [25j . The function / is a local free energy function taking into account 
short range interactions and entropy effects. The term involving J is the interaction energy due 
to a long-range, local mean field type interaction among the spins (particles), mediated by an 
interaction potential —^ d J(~fr) where 7 -1 is the range of the potential in units of the microscopic 
lattice spacing, and r is a lattice coordinate. This corresponds to the spin Hamiltonian 

H = ~l d J(l\x-y\)a(x)a(y) (1.3) 

x,y€Z<lf]j-^T L 

which specifies the long term attractive interaction. 
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The book [30] by Presutti may be consulted for more information on the relation between this 
functional and the underlying Ising system. Full details may be found there for the case in which 
the only short range interaction is a hard-core repulsion corresponding in the particle picture to 
the restriction that there is at most one particle at each lattice site. Then, with (3 denoting the 
inverse temperature, the local free energy f(m) is simply the —j3 times the lattice gas entropy term 

1 — m 1 — m 1 + m 1 + m 
s(m) = — log— — log— ^— ; (1.4) 

see [USD]. 

We shall focus on the case / = —(3~ 1 s here as well, though our analysis is readily adapted to 
the more general case of a strictly convex local free energy functional /(m); see [3] 

The mesoscopic length scale is set by the range of the J, i.e., 7 -1 lattice spacings. The macro- 
scopic length scale is L, the size of the domain. In the scaling limit that leads to .Fglp as a large 
deviations functional, the ratio between L and the range of J is held fixed (so that L is also mea- 
sured in units of 7 -1 , like the range of J), while 7 tends to zero. This is very different from the 
frequently encountered scalings in which one also takes a thermodynamic limit {L — > 00) at the 
same time one takes the continuum limit (7 — * 0), in which case the size L of the system grows 
faster rate than 7 _1 ; e.g., 7~ 2 . 

Our problem is to determine the minimum value of Fqlp subject to the constraint 
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m(x)dx = n , (1-5) 



and to characterize all of the profiles m that nearly minimize ^"glp • 

The nature of this minimization problem is clarified through use of the identity 

m{x) J(x — y)m(y)dxdy = - [ [ [m(x) — m(y)] 2 \ J(x — y)dxdy — f m 2 (x)dx , 



2 



which is valid for L > 1 (which we always take to be the case) on account of (|1.2p . to write the free 
energy functional in the form 



•F glp ijn) 



f{m(x)) - ^m 2 (x) 



dx + — I I [m(x) — m(y)] z J(x — y)dxdy . (1-6) 
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In case f(m)—m 2 /2 is a convex function of m, the minimization problem is trivial: One subtracts 
an appropriate multiple \im from f{m) so that the modified local free energy is minimized at m = n. 
Because of the constraint (|1.5p . this modification has no effect on the minimizers. Then, it is clear 
that the uniform profile m(x) = n is the unique minimizer. However, f(m) — m 2 /2 need not be 
convex. 

For example, with f(m) = —s[m)/j3, the function f(m) — m 2 /2 is strictly convex if f3 < 1, 
but for (3 > 1, it is not. Instead, it is a "double well" function with minima at imp, where mp is 
the positive solution to rap = tanh(/3m / g). In this case, ±7713 are the two equilibrium values of m, 
m±, mentioned earlier. To simplify the writing of T for this case on which we shall focus, let us 
introduce the function F on [—1,1] defined by 



F{m) 



1 < \ 1 2 
--s(m) - -m 
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this differs from —s{m)/f5 — m 2 /2 by a constant so that F(m) > 0, and F(m) = if and only if 
m = ±(3. 

We then define T to be 

T(m) := I F(m(x))dx + — f f [m(x) — m(y)] 2 J(x — y)dxdy . (1-7) 



The problem of determining inf^ <|jT(m) : J m{x)dx = n X and finding all minimizing 

profiles m is very easy for n £ [—1, — mp] and n £ [mp, 1]: In these intervals, the unique minimizer 
is the constant profile m(x) = n. (One easy way to see this is pointed out in Remark 14. 31 ) Hence, 
the question is only interesting for n £ (—mp,mp). 

This question has been well studied for fixed values of n £ (—mp, mp) as L tends to infinity (see 
[30| and references quoted therein). Here we are concerned with values of n such that 

d 

n = -771(3 + 0{L d + 1 ) , 

which turns out to be critical for droplet formation: If the value of n is low enough in this range, 
the minimizing profiles will be uniform, and if it is large enough in this range, the minimizing 
profile will represent a "droplet" of the +mp phase in a sea of the — mp phase; i.e., the droplet is 
the region in which m(x) ~ mp. 

In the well-studied problem, with fixed n £ {—7np,mp) and sufficiently large L, one always sees 
a droplet which has the equimolar volume Dq, specified by 

turDq - m (L d - D ) = nL d and hence D = - — V ^-L d . (1.8) 

2m/3 

Indeed, if a profile m takes on only the values imp then the constraint (|1.5|) will be satisfied if and 
only if the volume of the set on which it has the value +mp is Dq. 

d 

We shall see that when mp + n = 0(L <*+i) and L is large, then the droplets, when they exist, 
are smaller than this, but not too small: The volume D of the droplet will satisfy 

Dn < D < Dn . 



d+l 



That is, there is a universal lower bound on the size of stable droplets. Droplets that are "too 

2 

small" always "prefer to evaporate". The factor is independent of the particular interaction 

d+l 

d 

potential J, and if one considers n(L) defined by n(L) = —mp + KL d + 1 , one can determine the 
precise fraction of the equimolar volume as a function of K. 

Results of this type were first obtained for the 2-dimensional nearest neighbor Ising model by 
Biskup, Chayes and Kotecky [HIS]- In [5], they presented a general heuristic analysis of droplet 
formation, which predicts the above universal lower bound on D, and then they proved in [6] that 
the predictions of the heuristic analysis were correct for the nearest neighbor Ising model, by a very 
detailed analysis of the microscopic states that support the Gibbs measure. This rigorous analysis 
was carried out directly on the microscopic level, and did not involve the analysis of a free energy 
functional. However, their heuristic analysis, which was based on a competition between surface 
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tension and compressibility, could well be expected to apply to the GLP functional, and other free 
energy functionals like it. 

This expectation has been borne out in recent work of ourselves [11], and also of Belletini, Gelli, 
Luckhaus, Novaga [2]. In these papers, the problem of determining the nature of the minimizing 
profiles has been solved for a phenomenological analog of the GLP free energy functional, the 
Allen-Cahn free energy functional: 

1 I 6 2 I 

^(AC)(m) = - (1 -m 2 (x)) 2 dx + — / |Vm(a;)| 2 dx . (1.9) 

Here (1 — m 2 ) 2 /4 is a simple caricature of the more physical double well potential F(m), and the 
gradient integral penalizes variation in m just as 

\\ [ \m(x) - m(y)] 2 \J(x - y)dxdy (1.10) 
4 Jt l Jt l 

does. For purposes of this paper, we may regard Jyac) as a phenomenological caricature of GLP 
functional. Though J~(AC) is often called the Allen-Cahn, Cahn-Hilliard or Landau-Ginzburg 
functional, it goes back to van der Waals [T2j [32] . 

The analysis of F(ac) is much easier than the analysis of T because many tools that are 
applicable to the local interaction term J*^- | Vm(x)| 2 dx cannot be applied to its non-local relative 
(jl.lOp . In particular, both papers [IT] and [2] made use of an an idea originating with Modica [27] 
that used the co-area formula [J3] to get lower bounds on •7 7 (AC)( m )- The co-area formula is simply 
the change of variables in which one uses m itself as one of the variables of integration: 

d d x = 1 ., da h dh , 
| Vm{x)\ 

where do~h is the (d— l)-dimensional Huasdorff measure on the level surface Th = {x : m{x) = h}. 
Using this, 

^ c,(m) = L L (?' Vm(i)i + twM) dahdh ■ 

Then, by the arithmetic-geometric mean inequality, we have the lower bound 

JF (AC) (m) > f [ (~\l-h 2 \)da h dh 
> / e -{l-h 2 )\T h \dh, 

J[-l+e,l-e] 1 

(1.11) 

where e > 0, and |r^| is the (d — l)-dimensional Huassdorf measure of T^. Since for each h in the 
retained domain of integration, Th bounds a set including the set D e = {x : m(x) > 1 — e}, the 
isoperimetic inequality gives us a lower bound on \Th\ in terms of the ci-dimensional volume |-D e |- 
There is more to be done to determine the nature of the minimizers; see (1 1|, [2] . But the co-area 
formula, first used in this context by Modica and Mortola [28], provides a key to the analysis of 
the Allen-Cahn functional that does not seem to be adaptable to the non-local context where the 
interaction is not given in terms of a gradient. 
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On a technical level, our main innovations here are the development of tools to obtain a sharp 
lower bound in the non-local case. We shall make essential use of various rearrangement inequalities, 
combined with various truncation arguments and a priori estimates on near minimizers. The 
technical difficulties are worthwhile, we believe, because the GLP free energy functional has a direct 
connection with an underlying microscopic model, unlike the Allen-Cahn free energy functional. 

The rest of the paper is organized as follows. We first briefly recall the heuristic analysis in 
[5]. From this we deduce some natural conjectures about the nature of the minimizers of (jl.7l) . 
The conjectures are correct, so we state them as theorems, and prove them in the next sections. 
We conclude with a short section discussing further properties of the minimizers, and some open 
questions concerning them. 

2 The heuristic analysis 

Consider a free energy functional J-qlp as in (11 .6f) for some strictly convex function f[m) and some 
j3 > 0. Let g(m) = f(m) — (3m 2 /2, and suppose that this is a "double well potenetial" so that 
there is a tangent line am + b that touches the graph of g{m) exactly twice. Let (m_,g(m_)) and 
(m + , g(m+), m_ < m + , denote these two points. Then the functions g(m) — (am + b) is strictly 
positive except at m = m±, where it is zero. Notice that for any profile m(x) satisiying (|1.5|) . 
J Ti (am(x) + b)dx = L d (an + b), and so subtracting am + b from g(m) simply subtracts a fixed 
constant from 



and has no effect on the class of minimizers. 

Consider the case in which the constraint value n in (|2.ip is slightly above m_, and well below 
m + . There are several obvious options to consider when trying to construct profiles m(x) that will 
yield minimal, or nearly minimal, values in (|2.ip . 

One option might be to put all of the excess of m over m_ into a "droplet" in which m(x) = m+, 
and outside of which m(x) = m_. The constraint requires the volume of the droplet; i.e., the region 
in which m(x) ~ m + , would to be the equimolar volume Dq given by m + Do — m^(L d — Dq) = nL d . 

However, the optimal shape of the droplet would depend on the symmetry properties of J. In 
general, for L large compared to D^ d , one would expect the cost of forming a droplet to come 
from the surface tension between the m_ and m+ regions; see [31] , The shape that minimizes the 
surface tension is known as the Wulff shape for the functional ^"glp- Whatever the Wulff shape 
turns out to be, one would expect that for L large compared to D^ d , the free energy cost of forming 
a droplet of volume Dq should be proportional to Dq 1 ^ d . 

In particular, under our assumption that J is isotropic, we would expect the optimal droplet to 
be very nearly a ball of volume Dq, and the free energy cost of forming the droplet to be simply a 
multiple S, the surface tension, of the surface area of a ball of volume Dq. 

Therefore, letting denote the surface area of the unit sphere in M. d , we would have 




(2.1) 




(2.2) 
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for a profile that arranges the excess magnetization into a round droplet of volume Dq. Without 
the assumption that J is isotropic, the proportionality constant would be different, depending on 
the nature of the Wulff shape, but the cost of this option would still be proportional to Dq 1 . 

Another option might be to smear the excess uniformly over the background. This gives a bulk 
contribution, and its size is determined by the compressibility x— which is the inverse of </'(m_), 
where, as above, g(m) = f(m) — (3m 2 /2. Smearing a droplet of volume Dq over Tl, we get the 
uniform profile 

m(x) = n = rri- + (m+ — m— . (2-3) 
Hence, when L d is large compared to Dq, 

.AT" 



g{m(x)) « s(m_) + ^g"(m^ 



[mi — m_ , 



for all x. Integrating over the domain, we find that this option for m gives 

^GLP(m) — — d . 

Which option does better? The free energy in the droplet option is independent of L, while the 
cost of smearing the droplet over the backgound decreases as L increases. So with Dq held fixed, 
the droplet does better for small values of L, and the uniform profile does better for large values 
of L. To determine the break even point, equate the two values for T to find 

Sad ^ Dq _ ( m+ -m„) 2 Dl 

1- 



a d /dj 2x- L d " 



Thus, both are comparable when D^ d = 0(L d ), which by (|2.3[) means 



d 



n + m„ = 0{L d+r) . (2.4) 

This defines the critical scaling regime. 

What should one expect for the minimizing free energy in the critical scaling regime, and will 
the minimizers be given by some sort of droplet, or not? 

In [5], Biskup, Chayes and Kotecky proposed that to answer this question, one should introduce 
a volume fraction < rj < 1, and put t]Dq into a drop of the appropriate (Wulff) shape, and 
(1 — v)Dq into the uniform background. They then constructed a phenomenological thermodynamic 
free energy function <&(rj) which is the sum of the surface tension term and the uniform background 
term. In the case that J is isotropic, so that the Wulff shape is a ball, <&(n) is given by 

m = sJ^ lli + (m+ - m - ?({1 - n)D '' )2 



"did) 2 X - £<< 

The suggestion of [5] is that in great generality, one can resolve a competition between surface 
and bulk energy effects by choosing n £ [0, 1] to minimize <3?. This picture was then proven rigorously 
in [6] for the 2-dimensional Ising model with nearest neighbor interactions, where the Wulff shape 
is temperature dependent, and droplets correspond to connetcted clusters of + spins in a sea of — 
spins. 
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2.1 The reduced variational problem 

Define C(D ,L) by 



(m+-m_) 2 f d\ l/d Dl +1/d 



Then &{rf) can be written as 



D \ 1 - 1 / d 



For fixed Dq and L, the function <&(r]) is increasing at r/ = (with an infinite derivative), and 
has a local minimum for some 77 > 0. Depending on the value of C(Dq, L), the absolute minimum 
of $ on [0, 1] may be at either rj = 0, or at the local minimum at positive 77. In the first case, there 
is no droplet; everything gets smeared over the background. In the latter case, one has a droplet 
with a radius corresponding to the fraction of the excess that one puts in the droplet. 

To solve this minimization problem, notice that $(77) — ^(O) is a constant multiple of 

+ C(A>, L)V - 2C(L» , L)} , (2.6) 

which vanishes at 77 = 0, and hence has a minimum at some rj > if and only if it becomes negative 
somewhere. For which values of Dq and L does this happen? By the arithmetic-geometric mean 
inequality, 

> (^-^{(d + l)C V )^ (2.7) 
„ i d+1 

Thus, the quantity in (|2.6|) is minimized at 77 > if and only if 

C(D ,L)^^±± <2C(D ,L) . 
rfd+i 



Let C± be the value of C that gives equality in this last inequality. One finds, 

a 



d+l 

1 /d + l\~a~ 



d 

Moreover, with C(Dq,L) = C*, there is equality in the application made above of the arithmetic 

geometric mean inequality if and only if — ^ — r]~~3 = (d + l)C*r/. Therefore, define 77* by 77* = 

a 

d 

{dC+) d +! . One finds 

* = irn ■ (2 - 8) 

That is, if the minimum of <1> is attained at some positive value of 77, the positive value is never 
less than 77*. The volume of a droplet will always lie between that of the equimolar droplet Dq, 
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and the reduced value ti+Dq. Smaller droplets are never seen; they prefer to evaporate. Thus, the 
simple ansatz of dividing the excess volume between a droplet and the background, and optimizing 
over the volume fraction yields simple predictions for whether one sees a droplet of not, and the 
size of the droplet if there is one. 

Of course, the shape of the droplet was put into the ansatz by hand. The surface tension formula 

S ff d I — Ti I is what one expects to be appropriate for the isotropic long range interaction 
\Od/dJ 

that we consider here. For a non-isotropic interaction, one should have to replace this surface 
tension formula by the corresponding formula for the appropriate Wulff shape. The crucial point 
is that this would still be some multiple of Dq 1 ^ d , and would therefore lead in the same way to 
the analysis of rj[ri~~3 + C(Dq, L)n — 2C(Dq, L)], except for a different constant C(Dq, L). 

2.2 Statements of the theorems on the minimizers 

The theorems we present in this section refer to the specific GLP free energy functional T given in 
(jl.7p . This model has been particularly well studied in the literature (see |30| and the references 
there). As we have explained above, for n close to — nip the free energy of a round droplet of 
equimolar volume Dq should be given by (|2.2j) . and the ansatz of [5] suggests that 



inf iV(m) : ~T~rj [ m(x)dx = n\^ inf Sa d ( \r, l - l/d + C(D , L)(l - nf 
meM y L a J Tl J r?e[0.i] \cr d /dj I 

(2.9) 

where C(Dq, L) is given by (|2,5p using (jl.8p to express Dq in terms of n and L, and using the values 
of X- an d S appropriate to this particular form of T. Of course, X- can De easily computed from 
the explicit form of f(m) = — (l//?)s(m) — m 2 /2 given in (|1.4p . In fact, since s(m), and hence f(m) 
is symmetric in m, the compressibility is the same in the two phases, and we shall simply write x 
in place of X- for this model. 

As for the surface tension, an explicit variational formula for S is given in (|3.1|) in the next sec- 
tion. Suffice it to say here that it is the minimal value for a simpler variational problem concerning 
one dimension for profiles interpolating between — mp and +rng. 

We shall investigate (|2.9p in the critical scaling regime (|2.4p . for which n + rap is proportional 

d 

to L . Therefore, fix any K > 0, and define 



n = —mp + KL d + 1 . 



With this choice of n, (|1.8p gives us 



D = —L— and Sa d (—^-) = Sa d L~ . (2.10) 

2mp \<J d /dJ \2mpa d J 

Inserting this value of Dq in (|2.5p . one finds that C(Dq, L) = C(K), a constant depending only 
on K that is given by 

2ml ( d \ 3 / K \ 
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Therefore, with n = — rap + KL d + 1 , f|2.9[) becomes 



inf jjF(m) : — 7 / 



m(x)dx = n 



Kd V"^ 



L d+i 5^ inf 

2mp<J d ) r?G[0.1] 



r? i-i/rf + c(i ^ )(1 _ r?) 2 



Our first theorem says that if we divide both sides through by L d + 1 , this becomes exact in the 
limit as L tends to infinity: 

2.1 THEOREM. For all K > 0, 

d 2 -d ( If d 1 

lim L d + 1 inf < Tim) : — -, \ ra(x)dx = —ran + KL d + l > = 
L^oo m£M { L d Jr L J 

Kd 

where C(K) is given by \2.11\) . 



' v i-V* + C(K){l-ri) 2 , (2.12) 



By what has been explained in the previous section, the infimum on the right in (|2.12p occurs 
if and only if 



c(k) > c, 



d+1 

1 /d+l\~ 



d V 2 

Evidently from (|2.1ip . this is the case if and only if K > K* where 

d 

d+l ( xS\ d+1 fo-d\dTT / 01 „^ 

K * = — m \^) Vd) ■ (2 - 13) 

Theorem I2.ll therefore suggests that the curve 

n c (L) = -rap + K±L~^ 

is critical for droplet formation, so that for large L and densities n significantly below this level, 
the minimizers will be uniform, while for large L and densities n significantly above this level, the 
minimizers will correspond to spherical droplets of a reduced volume r/ c Do where rj c is given by 
(|2.8p . The following theorems bear this out. 



2.2 THEOREM. For all K < K„ and L sufficiently large, when 

d 

-1 < n < -rap + KL <*+i , 
the unique minimizer for T is the uniform order "parameter field ra(x) = n. 

To show that droplets do form for K > K*, we need a precise definition of we mean by a 

, d_ I 

"droplet of the +mp state in a sea of the — rap state". Toward this end, set k = [KL d + 1 J 3 and 
define the subsets A and C of Tl by slicing Tl at the following level curves of m: 

A = { x £ Tl : — rap + k < m(x) < rap — k } , C = { x 6 Tl : ra{x) > rap — k }. 
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2.3 THEOREM. Suppose K > K+, and n = —mp + KL~J+i . Then, given e > one can find 
a > so that for any trial function m for which T{m) < fi{ n ) + oc the following statements are 
true: 

||CUt1| - rj c D \ < eD 

and 

\\C\- Vc D \ <eD , 
where rj c is the optimal volume fraction from Theorem \2.1\ 

2.4 Remark. Since m(x) is close to, or larger than, mp on C , and is close to, or smaller than, 
—mp on Ti\{A U C), we may think of C as the "droplet of condensate" and A as the (evidently 
thin, and presumably anular) surface region of the droplet. The estimates of Theorem \2.3\ specify 
only the size of the droplet, and not its shape. The analysis of this shape is the subject of continuing 
research. It is heuristically clear from Theorem \2.1\ that this optimal shape is very close to a ball. 
Further discussion on what is required to prove this may be found in the final section. 

To prove the theorems we need good upper and lower bounds for .F(m) at admissible trial 
functions m. The upper bounds come from a trial function suggested by the ansatz in [5|. The 
lower bound is the part that is more technically challenging for the reasons explained above. We 
begin with the upper bound. 

3 The upper bound 

3.1 A good trial function for the intermediate regime 

The arguments of [5] suggest that one should use as a trial function a function of the form 

m^(x) = m {\x\ - r]dr ) + a(rj) , 

where: 

1. ro is the radius of a ball of volume Dq (known as the equimolar radius), and so r/dro is the 
radius of a ball with volume t/Dq. 

2. mo is some one dimensional transition profile that very nearly minimizes the cost in free 
energy of making the transition from m = —mp to +mp around the origin. 

3. a{rf) is a constant determined by the constraint / m/ ri \(x)dx = nL d . 

Do not confuse which is a functions on Tl, with mo, which is a function on M. 
• As 7] varies in the interval < n < 1, this family of "fractional droplet" trial functions interpo- 
lates between smearing everything over the background, for r\ = ; and putting everything into the 
equimolar droplet, for rj = 1. 



CCELM 



12 



3.2 Planar surface tension and the choice of m 

The natural choice for mo is given by minimizing the transition cost: 



S = inf 



1 



[f{m(z)) - f{mp)} + - / \m(z) - m(y)\ J(z~y)dy 



dz : lim m(z) 



where, writing x G M. d as x = (y, z) with y G R rf 1 and z G 



J(z) 



(3.1) 



(3.2) 



and J is the interaction potential in the GLP free energy functional. 

The quantity S is the planar surface tension; see |30| . It is well known [30] that the minimizer 
m is unique up to translations. In the rest of the paper, m is the minimizer vanishing at the origin. 

We now choose tjiq. We cannot simply choose mo = m since then mo(|x| — i] d ro) would not 
define a smooth, or even continuous, function on 7^. 

However, only mild modifications are required: We modify it so that mo(|x| — r/ 1//( Vo) defines a 
smooth function on 7/,, and the difference between mo and m goes to zero exponentially fast as L 
tends to infinity. We define mo (z) as any smooth function on R such that 



mo{z) 



d-l 

\z\ < Ld+1 



fh(z) if 

d-l 

-m ( gsgn(z) if \z\ > 2L d + 1 



d d-l 

By (|2.1U|) . in the critical scaling regime ro is proportional to L^ 1 , so that while L d + l is large, it is 
small compared to tq. 

3.3 The determination of a(r)) 

The constraint equation is / m^{x)dx = nL d , and since (by definition of ro) 

nL d = - mp [l° 

we have 



d r °) ' 



2m^r d -L d 



1 

77 



mo(|x| — r^dx . 



(3.3) 



We require sharp estimates on the integral on the right. But we know enough about mo to derive 
them, and can now estimate !F(m^ ) quite closely: 



3.1 LEMMA. In the critical scaling regime, with n = —mp + KL , 



n<n M )<L^r Sai{ —) 



n 



1-1/<Z 



+ C(K)(1- V f 



d —2d j ry 

Note that L d + 1 proportional to r ~ . In proving Lemma |3.1| it will be convenient to express 
our estimates in terms of powers of ro instead of powers of L. To prove Lemma 13.11 we start with: 
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d-l i 

3.2 LEMMA. For all rj such that L d +i < r rjd < r , 



/ mod 



x\ - r r]^)dx = mp 2^-r^-q - L d + 0{r d ~ 2 ) 



so i/mf 



a( ?? ) = 2m /3 ^^(l- ?? ) + 



Proof: Define = — mpsgn(z) and the constant M by M = J R z(p(z) — m(z))d,2, and set 
i 

r r\ = fo^S. Note that 



Tr. 



mo(|x| — r T] )dx = mp 



2 d r " " L . 



{p{\x\ - r v ) - m (\x\ - r v ))dx . 



Define 

h = 



\x\<2r v 



p(\x\ — r v ) — mo(\x\ — r v ))dx, and I2 



x >2r r 



- r,,) - m (|x| - r r) ))dx . 



We easily see that for all dimensions d, I2 = 0(e Ll/i ). Moreover, using polar coordinates 



h = <?d 



2r„ 



(p(s - r v ) - m (s - r v )) s d x ds 



Introducing the new variable z = s — r„, we see that if we extend the integration in z over the 
whole real line, we only make an error of size 0(e~ Ll 4 ) at most, and so 

d-l 



h = o: A r d - x Jjj>(z) - rn (z)) (l + dz + 0(e 



Taking into account the fact that (p(z) — mo(z)) is odd and rapidly decaying, we see that for d = 2 
or d = 3, 



d-l d—1 

(p(z) - m (z)) ( H ) dz = 



(p(z) — mo(z))zdz . 



' V / V 

In higher dimension this gives the leading order correction. This, together with the definition of 
7710(2) in terms of m(z), yields the bound on the integral. Then the bound on a{rj) follows from 
this and (13311. □ 



3.3 Remark. We see from Lemma \3. 2\ that in the critical scaling regime, except when i] = 1, 



a(rf) x U 



d+l 



(3.4) 



3.4 Computation of jF(m^)) 

With the trial function specified, we are now ready to prove Lemma 13.11 

Proof of Lemma I3.lt To simplify the notation, we write tuq to denote mo(|x| — r v ) and a to 
denote a(rj) so that m^) = tuq + a. We begin by estimating fj- F((m^)dx. Making a Taylor 
expansion, we find that for some A E [0, 1], 



F(m (v) ) = F(m ) + F'(m )a + -F"(m )a 2 + -F'" (m + Xa)a 3 . 
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We are required to produce a close upper bound on the integral of each of these terms over 
0,. It turns out that the terms with odd derivatives are negligible, and that to a very good of 
approximation 

F(muj\)dx ki I F(mo)dx H — J F"(mo)dxa 2 



f F(m )dx + — L d a 2 . 
JTr. 2 X 



(3.5) 

To see why this should be so, before going into the detailed calculations, note that F(—mp) = 

F(mp) = 0, and uiq is essentially equal to ±mg except in a shell of unit thickness and radius r^-qd . 

Thus, the term j r F(mo)dx is C(r^ _1 ). Likewise, since F'(—mp) = F'(jng) = 0, the integral of 

F'(mo) over Tl is ©(r^" 1 ). However, this integral gets multiplied by a, which is small. Hence this 

term is negligible compared to the first term. 

When we come to the second derivative term, we have F' (—mg) = F"(mg) = — and so the 

X 

L d 

integral of F"(mo) over is very close to — . Since this integral gets multiplied by a 2 , and we 

X 

have La x L >L^~ d ~+ I = O(r ~ ), this contribution it is of the same order as the first integral in 

the critical scaling regime. 

Likewise, we have an 0(L d ) bound on the integral of F"'(mo + Aq) over Tl, but while the 

integral involving F" gets multiplied by a 2 , this integral gets multiplied by a 3 , and so it too is 

negligible compared to the two integrals we shall keep. The next several paragraphs contain the 

precise calculations, and then we turn to the interaction term. 

1 1 -\- vci 

Note that F'(m) = 2^°^^ m an< ^ ~F'{±mfi) = 0- Since F'(m(z)) is an odd, rapidly 



decaying function of z, estimates just like the ones employed in the proof of Lemma (13.21) show that 
jf F'(m )dx = a d r d - x J F'(m{z)) (l + dz + Oie^") . 

Then, with the constant B defined by B = \ F 1 (m(z)) zdz , we have for d = 2 or d = 3 that 

Jr 



F'(m )dx = G d r d - 2 Brf^r + 0(e~ Ll/4 ) , (3.6) 



and in any dimension d > 4, the term <7d r o 2 Br] d gives the leading correction. 

Next, F"(mo) = -r- - 1 = — -, r. Therefore, with % = x(mg), 

(3(1 -m^) x("*o) 

/ \F"(m )dx = — —L d + R, \R\<(l-l)(l- v )2L,$. (3.7) 



Therefore the second order contribution is 

2ml 

-F"(m )dx = 



—a 2 / -F"(m )dx = — ^(1 - v)f + 0(L~W) 
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Finally, F"'(m) 



2m 



(1 — TO 2 ) 2 
1 



and so, for suffiently large L, 



'ft 1 



a 



—F"'(mo + Xa)dx 



< 



^L d a 3 = 0(L - *F 



°dr 



d-l 



(3.8) 



We now combine these estimates with an estimate on the interaction term. First note that by 
Lemma 13.21 the term — L d a 2 in (13. 5p is exactly the bulk term in <&o(rj). We shall combine the 
other term, fj- F(mo)dx, with the interaction terms, yielding: 

/ F(m )dx+j I I J(x - y)[m (x) - m (y)] 2 dxdy . 

It remains to extract the surface contribution to &o(rj) from these terms. 

We shall use the following simple fact: For any function g(y) depending only on \y\, and \x\ > 1, 
the range of J, we have 



J( x - y)g(\y\) d y 



l + O 



i 



i 



J(\x\ — s)g(s)ds . 



(3.9) 



A proof of such a statement in an even more general setting may be found in [17] . To prove the 
statement we need here, we simply change integration variables y \—* (\y\, z(y)), z(y) G R rf_1 in such 



a way that \y — x\ 



\y\ 



\x\) 2 + z(y) 2 . The change of variables doing this is explicitly given as 



follows: For s G M. + and z G R d 1 , set u = m and 



s ( u + |— r7(s, z) 



7/(8, z) 



(2\x\s) 



^l + l{s,z) 2 
Clearly \y\ = s and it is easy to check that \x — y\ 2 

l(s,z) = 



(2\x\s - \z\ 2 ) 2 



\ 2 . Moreover, 



uniformly in y in the unit ball around x, provided that \x\ is sufficiently large. 

Since J has unit range, we only need to bound the Jacobian of this transformation in the unit 
ball about x, and it is easy to see that for \x\ > 1 this Jacobian differs from unity by an amount 
that is uniformly bounded in the unit ball about a; by a multiple of (\x\ — l) -1 . 

Once more, estimates just like the ones employed in the proof of Lemma [3. 21 show that for some 
constant c, 



dx 



a d r d -\l + - 



F(m (x) + 7 / dyJ(x - y)[m {x) - m (y)} 2 
4 Jt l 



dz 



F{m(z))+ I dz'j(z-z)[m{z)-m{z')} 2 



where the errors are exponentially small in L 1 / 4 . But because mo is so close to m, this only differs 

2 

-J "1 1 1 1,1 si O s3 o ^* — 2o 

from Sad r o~ t] I by errors that are O(r ~ ). In the asymptotic scaling regime, r ~ x L d +! . 
Combining the estimates, we have the proof of the lemma □ 
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4 The lower bound 

The idea is, as in [10], to separate the surface and bulk contributions. The bulk estimate is similar 
to the one in [TO], while the surface estimate requires new ideas based on rearrangement arguments. 
The key to the lower bound is a partition ofl~L into three pieces: 

1. A region which will contribute a surface tension term to the free energy, 

2. A region which will contribute a compressibility term, 

3. A region that will make a negligible contribution. 

To do this we fix a number k > to be determined below. Define numbers h + and /i_ by 

= mp — k and /i_ = — mp + k . (4.1) 
Define the sets A, B and C by slicing Tl at the corresponding level curves: 



A 


= { 


X 




h- < 


m(x) 


<h+ } 


B 


= { 


X 


eT L 


m(x) 


< /i- 


} 


C 


= { 


X 




m(x) 


> h+ 


}■ 



(4.2) 

We denote by I a, Ib and Ic the contribution to J-{m) from the sets A, B and C. 
Define a radius R by 

^R d = \C\ , (4.3) 

where the right hand side denotes the measure of C. Evidently R is the radius of the ball with the 
same volume as C. 

It will be convenient in this section to write n in the form 

n = -mp + 5 . (4.4) 

d 

Notice that in the critical scaling regime, ixL d +i and from the definition of the equimolar radius 

5 = 2m^ = 2^| . (4.5) 
Given any trial function m(x), define fh(x) by truncating m(x) at the levels h- and /i+: 





V 


if 


??; 


» > K 


fh(x) = < 


m(x) 


if 


h- 


< m(x) < h 




h- 


if 


m 


» < h- 



(4.6) 



It is clear that 




\rh(x) — fh(y)\ J(x — y)dxdy < \m{x) — m(y)\ J(x — y)dxdy (4-7) 

' l ■ '' i J A J A 

< \m(x) — m(y)\ 2 J(x — y)dxdy . 
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Also clearly, 

F(m(x))dx = J F{m{x))dx + / F(m(x))dx + / F(m(x))dx . 

Jb Jc 



r L J a 



where F(m(x)) = f{m{x)) — f(mp). Therefore, 

F(m) > F s (m) + F B {m) (4. 



where 



4 



and 



!Fs(m) = - I I \fh{x) — fh(y)\ 2 J(x — y)dxdy + / F(m(x))dx (4-9) 
h L Jt l J a 

T B (m) = / F(m(x))dx . (4.10) 
Jb 

We shall refer to J-s( m ) as the surface contribution, and to ^^(m) as the bulk contribution. 
We shall obtain a lower bound on /l(w) by separately estimating these contributions. If the ansatz 
described in Section 2 is right, then at a minimizing m, essentially all of the contribution to T 
should come from Tsijn) and ^^(m), and so this lower bound will be quite sharp. 

In the next two subsections, we shall estimate Tsijn) and ^^(m) separately, starting with 
J-B(m)- First, however, we close this subsection by showing that if m is any trial function with 
J-{m) < J-(n), then C is not empty. (It is clear that if B is empty and k > 5, as we will assume 
later on, then the constraint J*^ m(x)dx = nL d cannot be satisfied.) 

For this purpose, it is advantageous to rewrite the free energy functional as follows: Define lu 
by uj(x) = m(x) — n. For m satisfying the constraint (|1.5j) . lu will satisfy 



Clearly, 



uj(x)dx = . (4.11) 



\m(x) — m(y)\ 2 J(x — y)dxdy = / \oj{x) — uj(y)\ 2 J{x — y)dxdy 
T L xT L JT L xT L 



Hence, if we define the functional Q by 



4 



where 



we have 



g(u>) = -[ \uj{x)-uj{y)\ 2 J{x-y)dxdy+ f G(u)dx , (4.12) 
't l xT l Jr L 



G(lu) = F(n + Lu) - F(n) - F'(n)cu (4.13) 



T{m) = T{n) + Q(u) , (4.14) 



since the term linear in lu drops out due to (14. lip whenever m satisfies the constraint (ll.5p . 
Thus, if J-{m) < J-{n), then Q{lu) < 0, which means that 



D := { x : G(lu(x)) < } / 



(4.15) 
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Figure 1: Plot of the function G(oj) 



The next thing to observe is that the set on which G(iv) < is a narrow interval (w-, w+) containing 
2mp whose width is of order 5 l l 2 . As long as we choose k large compared to 5 1 ^ 2 (we shall eventually 
choose k = <5 1//3 ), we will have D C C. 

The reason that this is true can be seen in Figure 1, where the function G(uj) is plotted. By 
()4.13p . one obtains G from F by subtracting the tangent line to the graph of F at m = n away 
from F, and then changing variables from m to to, which measures deviations from n. 

Since F is locally convex near m = n, G is locally convex near uj = 0, and so to = is one 
local minimum of G. Subtracting off the tangent line function "tilts" the graph of F downward to 
the right, and so the minimum of F at m = mp gets tilted to a slightly negative value; the global 
minimum of G lies in this dip below the axis. Since F'(n) = 0(5), the (negative) value of the global 
minimum is on the order of —2mpF'(n); i.e., also 0{5). But since the curvature at the minimum 
is strictly positive, the width of the interval on which G is negative is 0{5 1 ^ 2 ). 

The following lemma makes this, and more, precise: 



4.1 LEMMA. Let G be defined by \J^.13 ). Then the equation G(uj) = has exactly three solutions, 
0, CO- and uo + where < w_ < u+ . There is a constant c such that for all L sufficiently large 



\u± - 2m \ < c5 1/2 , 



(4.16) 



and, G(iv) < if and only if lo £ (w_,u;+). Consequently, if m is any trial function with ^(m) < 
J- in), and if k = S 1 ^ is used in the definition of and hence of C , then C is not empty. 

Moreover, G has a unique global minimizer uj+ which satisfies 



n 



+ oj+ - mp\ < c5 



for some fixed constant c depending only on F . If m is any trial function with ^{m) < T(n), there 
is an a < n so that 

"T-trunc(^) '■= max{ a , min{ n + w*, m(x) } } 

is also a valid trial function with .F(rra trunc ) < ^{m), and the sets A, B and C determined by mt rU nc 
are the same as those determined by m. 
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4.2 Remark. The last part of this lemma says that if we seek to prove any theorem concerning the 
sets A, B and C associated to a trial function m with J-(m) < T(n), then we may freely assume 
that m is bounded above by n + to* = mp + 0{5). It also implies that any minimizer m must be 
bounded above by mn + 0(5). 

Proof of Lemma I4.lt For the first part, we simply provide formulas that quantify the remarks 
in the paragraph preceding the statement of the lemma. It is clear from (|4.13p that uj = is one 
solution of G(u>) = 0. Also, since F'(n) = 0(5), and F(n) = 0(5 2 ), it follows that for some c < oo, 

G(u) > F(n + uj)-c5 . 

Thus, G(u) < requires F(n + uj) < c5; i.e., n + u must lie in one of the two "wells" of F. As 
before, let l/x denote F"(±mp), the second derivative of F at the bottom of the two wells. Let £ 
be defined by 

£ = iai \ m > : F"(m) > 



2x 

Then F(m) > F(t) > on [-£,£}, and so 

G{uj) > F(£) - c5 

for |n + u\ < £. For L large enough, F{£) — c5 > 0, and so we may restrict our attention to values 
of uj in the intervals — 1 < n + uj < — £ and £ < n + uj < 1. G is strictly convex in both of these 
intervals, and since G'(n) = 0, uj = is the unique minimizer of G in the left interval. 

To find the solutions in the right interval, introduce a new variable u defined u by n+uj =: —n—u 
and define H(u) := G(—2n — u). Then £ < n + uj < 1 if and only if — (n + 1) < u < —(n + £), and 
for such u 

H(u) = F(-n-u)- F(n) +2F'{n)n + F'(n)u 
= ^F"(-n-^)u 2 + 2F'(n)u + 2F'(n)n 

> —u 2 + 2F'(n)n + F'(n)u 

= — (u + A X F'(n)) 2 + 2F'(n)n-A X (F'(n)) 2 , (4.17) 

with £ £ [min{0, u}, max{0, u}] in the second line. We have used the Taylor expansion F(—n — u) = 
F(—n) — F'(—n)u + ^F"(— n — ^)u 2 and the fact that F(z) is even and F'(z) is odd, and finally, 
the lower bound on F" in the well. 

Evidently, H{u) > unless u_ < u < u+ where u± are the two roots of the quadratic expression 
on the right in (]4.17|) . Since F'(n) = 0(5), it is evident that there is a constant c such that 

\u±\ < cVo . 

By the local convexity of G, the remaining two solutions of G(uo) = must lie in the corre- 
sponding interval; i.e., (— 2n — it+, —2n — and G is positive outside this interval, except at 
uj = 
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When !F(m) < F{n), m is not constant, and / |w(x) — Lo(y)\ 2 J(x — y)dxdy > 0. Hence, 

JT l xT l 

for J-(m) < J~(n), we must have G(lo(x)) < on a set of positive measure. This proves the 
statements made in the first paragraph of the lemma. 

To prove the claims made in the second paragraph, note that G(lo) is an increasing function 
of lo to the right of its global minimum uj*. Since truncation always lowers the interaction energy 

/ \u(x) — Lo(y)\ 2 J(x — y)dxdy, replacing io(x) by min{ io(x) , lo* } lowers Q(uj). Note that, 
Jr L xT L 

if lo(x) was not already bounded above by a;*, the truncated function will no longer satisfy the 
constraint (|4.1ip . 

However, we can remedy this by another truncation at the other end: Let [lo]± denote the 
positive and negative parts of to. Then by (|4.11|) . 



tu}-dx = / [<jj] + dx 



and since 



a i — ^ / min{ \to\- , a }dx 
increases continuously from to /^[wj-dx as a increases from to 1, we can choose a so that 

/ min{ [lo]_ , a }dx = / min{ [lo] + , lo* }dx . 
Jt l ' Jt l 

Since G(lo) is a decreasing function of to on (— oo,0), this second truncation also lowers Q{lo), and 
restores the constraint (|4.1ip . Finally, truncating lo at —a and to*, corresponds to a truncation in 
m as in the lemma. The bound on \n + LO+ — mp\ comes from the fact that oo± is the unique non-zero 
solution to 

F'{n + Lo) = F'{n) , 

and what we have said in the proof of the first paragraph. □ 

4.3 Remark. As mentioned in the introduction, our minimization problem is trivial for n £ 
[—1,—mp] (and hence for n £ [7713,1] by symmetry). This is easily seen by considering the func- 
tional Q used in the previous proof: If n £ (—1,-771,3), then the tangent line being subtracted from 
the graph of F in \J\--13 ) would have a negative slope, and so subtracting it off would tilt the graph 



upward to the right, and not downward. Hence G will have a unique global minimum at lo = 0. 
Thus the unique minimizer of Q is the constant profile lo = 0, and then by \4- 14\ ) the unique mini- 
mizer of T is the constant profile m = n. For the values n = —1 and n = —mp, the situation is 
even more elementary. 

4.1 The bulk contribution 

The key to estimating Tb is that for k small enough, F is strictly convex on (— 1, /?,_): For any 

.,12 
h £ (—1, F (h) > — 1 + -jy- jp-- Define the quantity \- by 

— = f "^ = - 1 + -rt 1 w 

X- (3 I -hi 
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Then, by Taylor's Theorem, and using the fact that F{—mp) = F'(—mp) = 0, we have 

F(m(x)) > (m(x) +ma) 2 

2x- 



everywhere on {x € \ m{x) < hJ}. 
Therefore, 



F(m(x))dx 



> 



> 



B\ (j^J F(m(x))dx\ 

B\^~ f p J (m{x) + mp) 2 dx) (4.18) 
m(x)dx + mp\B\ 



1 

2X-IS 

This estimate should be quite sharp, since we expect any nearly minimizing profile m{x) to be 
nearly constant in B. Our problem is now reduced to that of estimating j B m(x)dx. Before going 
into the details, let us summarize what we would expect, and what lemmas we shall need, to prove 
what we would expect. 

First, we would expect the transition region A to be very "thin", so that: 

• \A\ is negligible compared to \B\ and \C\. 
In that case, we would expect 

m(x)dx ~ / m(x)dx — / m{x)dx 
b Jr L Jc 



Jc 



~ L d(_ mp + § )-^R* m p , (4.19) 

where in the last line, we have used the fact that m{x) is very close to m on C. In fact, on C, 
m(x) > mp — k by definition, and Lemma 14.11 and the remark following it will allow us to assume 
an upper bound of the form mp + 0{5). Also, if \A\ is negligibly small, 

\B\m « (L d - ^R d )m . (4.20) 
d 

Using (|4T20|) and ([TOD , we would have / m(x)dx + m \B\ ss L d (5 - 2m p ^(R/L) d ). Then 

Jb d 
using (|4.18p together with the simple (but not extravagant) bound \B\ < L d , we would have 

TB{m) * ^T d{5Ld ~ 2m ^ Rd)2 ■ (421) 

The next lemma gives the precise statement: 
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4.4 LEMMA. Let m be any trial function such that J-{m) < T{n), and such that m is bounded 
above by n + uj+, where lu± is defined and estimated in Lemma \4-1\ ■ Then 

« m) >£((^2^) 2 - £ ) , 

with e given by 

s = a\S- + f + c*)) (4.22) 



for some constant c. 

4.5 Remark. Note that if we choose k = <5 1/3 , and if R d /L d = 0{5), then e = 0(5 7 / 3 ). The 
surface contribution will limit the side of R, preventing cancelation in the main term, so that it will 
be 0(5 2 ), and hence strictly larger. 

To prove Lemma 14.41 we first need to show that \A\ is in fact negligible, as explained in the 
heuristics. The following Lemma takes care of that: 

4.6 LEMMA. Let m be any trial function such that J-{m) < J-{n). Then, for some finite c > 
and d > 

. A . Fin) 5 2 T H 
\A\<^-< c-L d . 

Proof: We have 

F(h+) = c'k 2 = F(h-) 

where d = ^F"(p), for some p with mp — k < p < mp and d > for k > small enough. It is easy 
to see from the definition of A and the properties of the function F that uniformly on A, 

F(m(x)) > F(h + ) = c'k 2 . 

Therefore 

I A > \A\c'k 2 . 

On the other hand, since J-(m) < J-(n), 

I A < F[n) = F(n)L d . 

Since F(n) = c\5 2 , we get the result. □ 

Now that we have Lemma 14.61 we return to the proof of Lemma 14.41 
Proof of Lemma 14. 4t Note that 

™ W dx = ^-/» W dx-/ m (x)dx. 

B Jc J A 
By Lemma 14.61 

-mp\A\ < h_\A\ < [ m(x)dx < h+\A\ < mp\A\ 

Thus, 



.4 
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/.nMd.-M/'-ICHO <m, W+ | C |m,-/m(,)d 
Since n = — m« + (5 and it is evident that \B\ = L d — ^R d — \A\, 



[J m(x)dx + mp\B\\ - (6L d - 2mp\C\) 



d 

< 2m /3 \A\ + 



|C|rng — / m{x)dx 



c 



(4.23) 



Next, on C, m(x) > mp — k by the definition of C. Also, by the hypothesis that m{x) < n + o;* 
for all x, and Lemma 14. lj, m(x) < mp + c5 for some fixed constant c, and for all x. Thus, 



|C|rn(3 — / m{x)dx 
c 



< \C\(k + c8) 



and hence, (14. 23ft yields 

J m(x)dx + mp\B\^J - (5L d - 2m p \C\) 



< 2mp\A\ + \C\{k + c5) . 



Since the inequality \a — 6| < c implies a 2 > b 2 — 2\b\c, we have 

J m(x)dx + m^l-BM > 

(<5L d -2m /3 |C7|) 2 -2 5L d -2m /3 |C| (2m/}\A\ + \C\(k + cS)) . (4.24) 
Going back to (|4.18j) and using the estimate \B\ < L d and the definition \C\ = adR d /d, we obtain 

L d ( . . <r d fl'" ' 



2X- 



5 — 2mp- 



d L d 



\ A \ VdR d , 
2m(3 T d+ ~dl? {K+ } 



Now using Lemma 14.61 to estimate \A\, we obtain the result. 
4.2 The surface contribution 

Our goal in this subsection is to prove the following estimate. 

4.7 LEMMA. Let m be any trial function such that J-(m) < J-(n) 

1 



(4.25) 
□ 



1-0 



R-2 



a d [R-2] 



d-l 



1 



rap) 



(4.26) 



where S is the surface tension, and [a]+ = max{ a , }, so that the bound is trivially true for 
R < 2, twice the range of J. 
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Let us first explain the heuristics, and then collect the lemmas required to substantiate them. 
To prove this lemma we need to relate J-g to the one dimensional functional defined in (|3.ip . which 
gives the planar surface tension S. To do this, we use rearrangement inequalities to replace our 
near minimizer m by a radial function on all of M 2 . This radial function will give us a trial function 
for M . 

The Riesz rearrangement inequality that we intend to use applies to functions on W d , and not 
on the torus, hence the first thing we have to do is to extend J-s to a functional on profiles in all 
of M. d without lowering the value of J-s too much. Here is why we can expect that this is possible. 

We expect that for a non-constant near minimizer m, C should be essentially a sphere of radius 
R that we can take to be centered in Tl, considered as a <i-cube of side length L in W 1 , and that 
the whole transition region A will be in an annulus close to C. In particular, the truncation m of 
m that is defiend in (|4.6p and used in the definition (|4.9j) of J-s satisfies fh(x) = h- for all x within 
unit distance of the boundary of the square. Now extend in to a function m on all of M. d . Do this 
by defining m(x) = h- for x outside Tl; i.e., 



mix) 



fh(x) if x G Tl 



if x € 



(4.27) 



L ■ 



Then, since J is supported by the unit sphere, it would follow from m(x) = /i_ everywhere near 
the boundary of Tl that 



\fh(x) — rh(y)\ 2 J(x — y)dxdy 
As for the potential term, define F by 



\m{x) — m(y)\ 2 J(x — y)dxdy . 



(4.28) 



F(m) 



F(m) if h + > m > /i_ 

if m<L ormXiJ 



(4.29) 



With this definition, 



F(m(x))dx = / F{fh(x))dx 
A J A 



F(m(x))dx . 



(4.30) 



Then combing (gSD, (|4T2"g|) and (|4T30]) . we would have 

1 



J^si 111 ) > / F(m{x))dx + 



m(x) — m{y)\ J(x — y)dxdy . 



We are now in a position to use rearrangement inequalities to make contact with the one 
dimensional variational problem f)3. 1 j) that defines the planar surface tension: Let m* denote the 
spherical decreasing rearrangement of fh (see [23]). Then by the Riesz rearrangement inequality 



\m(x) — m(y)\ 2 J(x — y)dxdy > 



\m*(x) — m*(y)\ 2 J(x — y)dxdy , 



and of course 



F(m(x))dx = / F(m*(x))dx 
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Therefore, if our intuition about the size and shape of C is right, we should have 

1 



T s {m) > 



\m*(x) — m*(y)\ J{x — y)dxdy + / F(m*(x))dx . 



Let r = \x\. Because of the spherical rearrangement, m*(x) depends on \x\, and the corresponding 
region C is indeed a sphere of radius R. 

The key to making contact with the planar surface tension is the fact that m*(r)/(l — ^-), 
extended by mp for r < 0, is a valid trial function for the one dimensional variational problem (|3.ip 
defining S. We shall use this fact to get a lower bound on !Fs(rn) for a non-constant minimizer 
m that is of the form (|4.26p . Note that apart from some small corrections, the main term in this 
bound is a c [R d ~ 1 S, the contribution we would expect for a droplet of radius R. 

To carry out this program of estimation, we need to show that \C\ is not too large: If \C\ is 
large, then it is easy for C to "wrap around" so that (I4,28h is not even approximately true. Then 
we would be prevented from applying the Riesz rearrangement inequality. Then next lemma shows 
that in fact, if m is a trial function with !F(m) < J-(n), then C, which we know to be non-empty 
by Lemma [4TT| has volume \C\ < 0(L d 5). 

4.8 LEMMA. Letm be any trial function such that F(m) < F(n). Then, for L sufficiently large, 
there is a constant c such that 

|C| < cL d 5 = cL^ . 

Proof: Since 

nL d = (—mp + 5)L d = / m(x)dx + / m(x)dx + / m(x)dx, (4-31) 

J A Jb Jc 

we can use the obvious lower bounds m(x) > —1 on A U B and m{x) > h + on C to conclude that 

(-mp + 8)L d > -(\A\ + \B\) + (mp - k)\C\. 
By using \C\ = L d - (\A\ + \B\), we get 

(|^| + |B|)(1 + m p - k) > {2mp - 5)L d . 
By Lemma 14.61 \A\ < c5 2 K~ 2 L d . Hence, for L sufficiently large 

\B\ > -^L d . 

1 + 771/3 



On the other hand, by (|4.18p 

Tin) > !F(m) > [ F(m)dx > - — - ( j=t [ m(x)dx + ran ] . 

Jb 2 X- \\B\ Jb J 



Therefore 



1 



\B\ JB 
by using !F{ri) < c±S 2 L 



m(x)dx + 777/3 



B 
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w 



Figure 2: Possible droplet shape. 



Finally, using this in (|4.3ip we have: 
(-m p + 8)L d >-(L d -\A\ 

Then 



\c\y 



B\ 



m(x)dx — \A\ + (rri/3 — k)\C\. 



\C\(2mp — k — c 2 5) < L d (5 + c 2 5) + \A\. 

Using Lemma 14.61 to bound \A\, and taking L sufficiently large we conclude the proof. □ 
Armed with this lemma on \C\ , we return to the proof of Lemma [4.71 In our heuristic discussion, 
we relied on our expectation that C is nearly a disk centered in Tl (with an appropriate choice 
of the origin in Tl) with a radius small compared to L in order to justify (|4.28p . At this stage 
something less can be proved, which still suffices for the proof of Theorem l2.lt 

•If m is a trial function with J-{m) < J-(n), and if C does "wrap around the torus" Tl, then the 
arms that "wrap around" are very thin, as shown in Figure 2. 

For any of the d coordinate directions 1 < i < d, consider the volume of A U C that is contained 
in the slab a < x\ < a + 2. If for each choice of a this volume is at least w, then 

\AUC\>^ W , 

since integrating this volume in a from —L/2 to L/2 gives twice the volume of \AL) C\. Hence there 
is at least one choice of a for which 

2|AUC| 



\{x : a < Xi < a + 2 } D (AUC)| < 



L 



Now by the translation invariance of we may freely translate m, and so may assume that a 
Thus without loss of generality, we may assume that for each coordinate direction i, 

2\AUC\ 



{ x : \xi±L/2\ < 1 }n(AUC) \< 



L 



L/2. 



(4.32) 
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4.9 LEMMA. There is a constant c so that for any trial function m with !F(m) < J-(n), 

I I \rh(x) — fh{y)\ 2 J(x — y)dxdy > f f \m(x) — rh{y)\ 2 J(x — y)dxdy — 2d} 
and for some finite c, 



L 



JAUCI 
2d- — - — - < cL 



Proof: The integral on the left hand side can only be smaller than the integral on the right hand 
side only on account of pairs of points (x, y) with, say, x in the cube representing 7^, and y outside 
it, and where x £ A U C, since otherwise m(x) = m(y) = h-. Since J has unit range, x must have 
\xi rb L/2\ < 1 for at least one 1 < i < d. Hence the total contribution from such pairs of points is 
a fixed multiple of \A U C\/L, by (|4.32p . Then using our bounds on \A\ and C from Lemmas 14.61 
and 14.81 we obtain final bound. □ 



Lemma [4.91 gives us the rigourous replacement for (|4.28p in our heuristic analysis. We now turn 
to the term involving F. Note that \x\ > R at all points x in the support of F(m*(x)). Hence, 
going to spherical coordinates, 



F(m*(x))dx = a d / F{m*{r))r d - l dr > a d R d - L F(m*(r))dr 
Jr Jr 

To proceed, we once again use (j3.9|) . By construction, m*(r) = mp — k for all 

i 

d 



(4.33) 



and the r.h.s. is larger than R. 
Therefore, if r < R — 1, then 



for all s, since J has unit range. 
Hence by (l33|) . 



r I -(\A\ + \C\) 



\m*(r) — m*(s)\ 2 J(r — s) = 



(4.34) 



\m*(x) — m*(y)\ 2 J(x — y)dxdy 



1 - O 



1 



R-2 



\m*(r) - m*(s)\ 2 J{r - s)a d r d ' l drds , (4.35) 

where once again, [•]+ is the positive part function. Indeed, the left hand side is clearly positive, so 
if R is sufficiently small that using the uniform bound on the Jacobian that led to (13. 9|) provides a 
negative lower bound, simply use the trivial lower bound by zero instead. 

Next, since (|4.34[) implies J(r — s)r d_1 > [R — l] d T 1 J(r — s) for all r, s > 0, we have 



\m*(x) — m*(y)\ 2 J(x — y)dxdy 



1 - O 



1 



R-2 



\R-l 



d-l 



oo roo 



JO 



m*(r) — m*(s)\ 2 J(r — s)a d drds , (4.36) 
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We are now ready to prove Lemma l4.7t 

Proof of Lemma 14. 7t By (j4.30p and (|4.33[) for the terms involving F, and Lemma l4*?9l and (|4.36p . 
it remains only to show that 

- / \m*(r) -m*(s)\ 2 J(r - s)drds + / F(m*(r))dr > (1 ) 2 S . 

4io io Jo mp 

Our proof of this rests on the fact that m*(r)/(l — ^-), extended by mp for r < 0, is a valid trial 
function for the variational problem defining S. To show this, we need to prove that it approaches 
the correct asymptotic values; i.e., that it satisfies the constraint imposed on (|3.ip . In fact, due to 
the lemma UTTl C in not empty when T{m) < J-(n), and of course B is not empty too. So, by the 
definition of C, m*(r) = mp — k for all r < R. Consequently, m*(r) — m*(s) = if r, s < R so that 
provided R > 1, the range of J we can extend the region of integration to (— oo, +oo): 

\m*(r) — m*(s)\ 2 J(r — s)drds = 



o 

(1-— ) 2 / l\m*{r)/(l-—)-m*{s)/{l-—)\ 2 J{r-s)drds. (4.37) 
mp J R J R mp mp 

Now, if R < 2, the bound in the Lemma is trivial, and there is nothing to prove, Hence there is no 
harm in assuming that R > 2, which we now do. 

Since F is decreasing on m > 0, and since for m > 0, m/(l — ^-) > m, it follows that 



F(m*(x)) > F m*(x)/(l ) > (1 ) J F m*(x)/(l ) 

m/3 / mp V m /3 



Therefore, 



oo />oo 

2~f 



JO 



m*(r) — m*(s)\ J(r — s)drds + / i ? (m*(r))dr > 



7(!-— ) 2 



K 



\m*(r)/(l ) - m*(s)/(l )|V(r - s)drds (4.38) 

mp mp 



+ / F(m*(r)/{1- —))dr 



> ( i__!L)25. 

The proof of Lemma 14.71 is complete. □ 



5 Proofs of the theorems. 
5.1 Proof of the Theorem 12.11 

We fix the value n = 5 1//3 for the proof. Given a trial function m with J-(m) < J~(n), we replace 
m by its truncation as defined in Lemma 14.11 This lowers the free energy, and does not change 
the sets A, B and C. In summary, after this replacement we have a trial function that has a free 
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energy at least as low as the one we started with, the same value of R, and which is bounded above 
by n + <jj+ as in Lemma 14.11 

Then using Lemmas 14.71 and 14.41 we conclude that 



> (yd 



1-0 



l 



R 



[R - 2]*" 1 (1 - S + ^-(6 - 2m^(R/L) d ) 2 - ^-e . 

(5.1) 



It now remains to optimize this over R. We note that the lower bound on J-B(m) decreases as R 
increases until R is of order tq. But for such values of r, the expression for the lower bound on the 
surface contribution simplifies: 

1 



1 - O 



R-2 



[R-2]^ 1 = R^ + Oi^- 2 ) 



Now introduce S- = S[ 1 ) and ri = R d lr d . Then we can rewrite this lower bound as 

Hm) > S-ajrt 1 + j^C(n)(l - ^ - ^-e + 0{r d ~ 2 ) . (5.2) 

By LemmalMl \A\/L d = 0(5 2 /k 2 ), and by LemmalS R = O(L^r). Thus, (R/L) d = 0(8). 
Therefore, 

£ = O ( ~ + 5 



With the choice k = 5 1//3 , this gives us e = 0(5 7 ^ 3 ), as noted in the remark following Lemma 14.41 
The essential point is that this is negligible compared to 5 2 as L tends to infinity in the critical 
scaling regime. 

As L — ► oo in the critical scaling regime, S_ — > S and x~ ~^ X- Moreover, from (|4.5p and 
the definition of e, L "e/cr^r _ = L 3 ( d + 1 ) — > as L — > oo. Thus, (15. 2D provides the lower bound 
needed to prove (|2.12|) . The upper bound is provided by Lemma [37TJ The remaining statements 
follow from the analysis of the minimization of the phenomenological free energy function &(rj) that 
was explained in Section 2 □ 

5.2 Proof of Theorem [2721 

Suppose that n = —mp + KL d ^ d+1 ^ where K < K+. We shall show that, in this case, any non- 
constant trial function m has a higher free energy than the uniform trial function m(x) = n, at 
least for all sufficiently large L. 

Recalling that Sad,r d ~ 1 C(Do, L) = !F(n), define fj by 

77 = su P |r ? : S^" 1 L 1 ' 1 ^ + -J^-C(D , L)(l - V A - ^-e < S<r d r d - l C(D ,L) \ 

As in the proof of Theorem 12.11 for all L sufficiently large, S- is sufficiently close to S, and x— 
is sufficiently close to x that 
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for some C < C*. For C < C*, the unique minimizer of 

77 i-c 7/ 1_1/d + C(l - V? 

is 77 = 0. Therefore, since eL^/cr^o" 1 - * as L — > 00, it follows that 77 — > as L — > 00. 

Now, as in the previous section, for any non-uniform minimizer m, there is a relation between 77 
and the size of the level set \{m > mg — «}|: for a given 7? = (R/ro) l / d , \{m > mg — k}\ = (ad/d)R d . 
Here, as in the last section, k = with 5 given by (|4.5I) , It follows from (15, 2j) and the definition 
of 77 that for any non-constant minimizer m, 77 < 77, and so \{m > mg — k}\ is negligibly small 
compared with Dq = (<Jd/d)rQ, the volume of the equimolar ball, when L is large. 

In other words, if n = —mg + KL d ^ d+1 ^ where K < K*, and L is large, then any droplet in any 
minimizer must be extremely small. To prove Theorem 12.21 it therefore suffices to show that such 
extremely small drops are impossible in a minimizing m. The following Lemma gives the required 
lower bound, and completes the proof. □ 

5.1 LEMMA. For all K > 0, there is a constant ck > depending only on K so that if n < 

d 

— mg + KL d + 1 and m is any non-uniform minimizer for then 

^j-R d := \{x e T L I m(x) >m 3 - k}\ > c K r^ . 
Moreover, ck is uniformly strictly positive for all K in an interval around K+. 

We first explain the idea behind the proof. It is convenient to use (|4.14p to write J-(m) in terms 
of G{oo) as in the proof of Lemma 14.11 

We know from Lemma [4. 11 using the notation from there, that if !F(rn) < J-(n), then Q{oS) < 0. 
Since for some c < 00, G{ui{x)) > —cS for all x, and since the set on which G(u>(x)) < is contained 
in C, 

dy\oj(x) — oo{y)\ 2 J(x — y)dxdy + / G(ui(x))dx 



G{uj) > 



4 



A J A 



c5R d . (5.3) 



Now, if we could show that for some S > 0, 

1 

4 




\u(x) — u(y)\ J(x — y)dxdy + / G(co(x))dx 
IAJA J A 



we would have 

g(u) > ^j-SR' 1 - 1 - c5R d = R^ 1 {^S - c6R) , 

and this is strictly positive unless R is on the order of and in the critical scaling regime, ro is 
proportional to 5~ 1 . 
Of course 



\uj(x) — uj(y)\ 2 J(x — y)dxdy + / G(u}(x))dx 

! \ ■/.! J A 

1 



4 



m(x) — m(y)\ J{x — y)dxdy+ I G(m(x) — n)dx . (5.4) 
A J A J A 



so we can hope to bring the methods of Section 4.3 to bear on this problem. There are two obstacles: 
The first obstacle is that the possible penalty for "unwrapping the torus" that is estimated in 
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Lemma HT9l is of order (\A\ + \C\)/L. When R is small, so is \C\, but the only upper bound that 
we have on \A\ is the one provided by Lemma 14.61 This is independent of R, and so it very well 
can be that for small R, \A\/L is large compared to so that the penalty for "unwrapping the 

torus" completely swallows up the surface term. 

The second obstacle is that lower bound on the surface contribution that we obtained in 
Lemma 14.71 becomes trivial for R < 2, twice the range of J. 

To deal with the first obstacle, use the fact that on A, G(m(x) — n) > ck 2 for some c > 0. 
Therefore, 

/ G(m(x) — n)dx = — / G(m{x) — n)dx H — |^4|ck 2 . 
J A 2 J A 2 

We then have from (|5.4|) that 




A J A 



Since 



\uj{x) — uj{y)\ 2 J{x — y)dxdy + l G{uj{x))dx = 

J A 

- \m(x) — m(y)\ 2 J(x — y)dxdy + - / G(m{x) — n)dx + -| A\ck 2 . (5.5) 
4 JaJa 2 J A 2 

2 \A\ck > L 

for all large L, the final term in (15. 5ft more than compensates for the price of "unwrapping the 
torus" . 

It remains to deal with the second obstacle, as concerns 

- / [ \m(x) — m(y)\ 2 J(x — y)dxdy + - f G(m(x) — n)dx . 
4 J A J A 2 J A 

Toward this end, we first give a simple, direct argument, independent of the reasoning in Section 
4.3 to show that if ^(m) < Tin), and m is not constant, then R is bounded below by a constant 
of order one. Specifically, we shall show that in this case , R > 2~ <yl+l l d \ 

We then "cut down" the range of J so this it is small compared with 2~^ 1+l ^ d \ That is, for 
p > 0, define J p (s) = J(s) for < s < p, and J p (s) = otherwise. Having reduced the range of 

J, the error term in (|3.9p becomes ( 1 — O ( — — — ) ) for |x| > R — 2p. We choose r sufficiently 

V \ r - 2 pJ J 

small that this factor is at least 2/3 for R > 2~^ 1+1 / d \ This shall take care of the second obstacle. 
We now provide the detials. 

Proof of Lemma I5.lt We shall first show that if J-{m) < J-{n), and m is not constant, then 
R > 2 - ( 1+1 / d ). First of all, by Lemma |4. H if !F{m) < J-(n), and m is not constant, then C is not 
empty, and at least we know R > 0. Let us improve on this. 

With uj_ defined as in Lemma 14. 1\ define C := { x : u>(x) > u>- }, and define R so that 

^-R d = \C\. Then by Lemma HH there is a constant ci so that / G{uj{x))dx > —c\5R . 
a Jt l 

We now claim that if R < 2~( l+l l d \ the interaction term makes a much larger positive contri- 
bution to G(u>), so that G(u>) > 0, which would imply ^{m) > F{n). 

To see this, note that by Lemma |4.1| if x E C, and y ^ C, then for some positive constant c, 

\m(x) — m(y)\ 2 = \lu(x) — u>(y)\ 2 > en 2 = c5 2 ^ , 
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since m(y) < mg — k, but m{x) > mp — 0(5 1 ^ 2 ). 

Now, by our assumptions on J, for each x S C, J(x — y) > a > on the ball of radius 1/2 
about x. If R < 2-( 1+1 / rf ), then C can fill up at most one half of this ball, and so the volume of 
the set of points y in the ball of radius 1/2 about x for which y ^ C is at least < ^j-2~ <yd+l \ Hence 



da; 



a 



dy\uj(x) - uj{y)\ 2 > Idx dy\w(x) - u(y)\ 2 J{x -y)> c 2 S 2/3 R° 



for some c 2 > 0. Altogether, Q(u) > [c2<5 2//3 — c%S]R d , and this is strictly positive for L large enough. 
Therefore, we may assume without loss of generality that m is such that R > 2~^ l+l l d \ 

We return to (|5.5p . and "cut down" the range of J, replacing J by J p , p < 1, which clearly 
decreases the right hand side of (|5.5p 




|m(x) — m(?/)| 2 J(x — y)dxdy + 



A J A 



If 1 

- / G(m(x) - n)dx + -\A\ck 2 > 
2 J A 2 



7 / / |?n(a;) — m(y)| 2 J p (x — y)dxdy + - / G(m{x) — n)dx + -| A\ck 2 . (5.6) 

^ JaJa % Ja 2 



Next, by (|4.7p and Lemma 14.91 we have 
m(a;) — m(y)| 2 J p (x — y)dxdy > 



A J A 



\m(x) — m(y)\ 2 J p (x — y)dxdy — 2d 



\A\ + \C\ 
L 



■ (5-7) 



Now let J p be defined in terms of J p just as J is defined in terms of J, and let m* denote the 
the rearrangement of in, as before. Having reduced the range of J from 1 to p, the estimate (I4.36P 
becomes 



\m*(x) — m*(y)\ 2 J p (x — y)dxdy 



1 - O 



P 



[R - p]_ 



m*(r) — m*(s)\ 2 J p (r — s)a<idrds , (5i 



R-2p 

where once again, [•]+ is the positive part function. 

We now choose < p < 2~( 2+1 / rf ) small enough that the prefactor 



is 



at least -R d 1 for R > 2 Of course, we also have 

2 



1 - O 



P 



R-2p 



G{u{x))dx > OdR 11 - 1 / G{m*{z) - n)dz 
where as in ()4.29p . we define G by 

\G(u)) if h+ > lo + n> h- 
G(u>) = < 

10 if w + n<L orm>/i + 

Combining (joUj) . ([5T7|) . ([STHjl and ([5\9]) . we have 



[R-pl 



(5.9) 



(5.10) 




1 



m(x) — m(y)\ J(x — y)dxdy -\ — / G(m{x) — n)dx H — |^4|ck > 



A J A 



1 



m(r) — m(r)| J p (x — y)addrds + — / G(m(r) — n)dr 



\A\ + \C\ 1, 



+ -|^|ck 2 

2L T 1 



(5.11) 
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(5.12) 



Now let S be defined by 

S = inf < — / dz [ G(m*(z) — n) H — / \m(z) — m(z')\ 2 J p (z — z')dz' ] : lim m(z) = 
[ 2 J R V 4 / z-»±°° 

Combining (|5.3p (|5.5p and (|5. 1 If) . we finally have 

g(o;)>5^-d^ + ( ; r ^ + icK 2 |A|. 

n 2d Ml 1 

Since k 2 = L 3(d+i) ; f or i arge ^ -2-|— 1 + -ck j |^| > 0, and 

<7( w ) > s^-i - - ctfif = R d ^ (s^-c^- c5R 

a dL \ a dL 

This is positive unless R > c5 _1 , and by (|4.14p we conclude that whenever J-(m) < J-(n), and m is 
not constant, R > cd^ 1 . Since in the critical scaling regime, ro is proportional to this proves 
the lemma. □ 



5.3 Proof of Theorem IX3l 

Suppose m is such that ^{m) < /x,(n) + a, for some given a > 0. With the notation of the previous 
sections, let rj = \C\/Dq and r\ c be the optimal volume fraction corresponding to n. Then, by (15. 2p 
we get 

F{m) > S- [$(r?) -$(r? c )] + 5_$(?7 C ) + o(LSTT) (5.13) 
> 5_[*(77)-*(77 c )]+/x(n) + o(L3TT). 

d 

This entails that $(r/) — $(7? c ) < #- + o(L d +! ). By the definition of $, there is a constant </?o such 

|C| 



that ^o(?7 — Vc) < o-^ d+1 + o(l), so, for L sufficiently large 



6 The shape problem 



n 11 



<o(l). □ 



Throughout this section, we use the notation defined in (|4.ip and (|4,2p . In particular, given a trial 
function m, the set C is the set of points x on which m(x) takes values that are close to mp or 
larger, B is the set of points x on which m(x) takes values that are close to — or smaller, and A 
is everything else. We shall also assume that m is defined on all of we have already seen how 
to extend m from Tl to M. d with negligible cost in free energy, so let us suppose this is done, and 
\A\ and |C| are finite. Finally, in this subsection we use the notation (|4.3p . so that R is the radius 
of the ball in M. d with the same Lebesgue measure as C. Up to now, we have been concerned with 
the sizes of A and C. Going forward, we are concerned with their shape. It is easier to get control 
of this in the local case. 

In the case of the local (Allen-Cahn or van der Waals) free energy functional (jl.9p . the lower 
bound (jl.lip brings the surface area of the boundary of C into the lower bound on the free energy. 
Then stability results for the isoperimetric inequality [T9| [20] can be used to show [Tl] that if C 
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is not nearly spherical, there is a significant cost in free energy. It is an open problem, which we 
refer to as the shape problem, to prove this for the Gates-Lebowitz-Penrose free energy functional. 
To clarify the difference between the local and non-local cases, let us briefly recall an argument 
from [llj. (We shall in fact improve the result in [11] by using a new stability inequality from [16J.) 

Let E be a Borel measurable set in M. d . As usual, let \E\ denote is Lebesgue measure, and 
let P(E) denote its perimeter. If the set E is sufficiently regular, P(E) is the d — 1 dimensional 
Haussdorf measure of the boundary of E, though for very irregular sets it can be much smaller. 
The perimeter functional is an extension of the surface area functional to general Borel sets, due 
to De Giorgi [13J , with the key property that it enjoys good lower semicontinuity properties. See 
Maggi's review [26] for more information. 

The isoperimetric deficit of E, 6(E), is the quantity 

8(E) = ^ 1 . (6.1) 

d {d-i)/d a Vd\ E \{d-\)/d 

The general isoperimetric inequality of De Giorgi says that 5(E) > with equality if and only if 
E, up to a set of measure zero, is a ball. 

The stability results that we refer to give a lower bound on 6(E) in terms of the Fraenkel 
asymmetry of E, A(E), which measures the extent to which E differs from being a ball: 

A(E) = inf { |gA g r>iC)l : = \E\ , x G R d } (6.2) 

where B(r, x) is the ball of radius r centered on x in and EAB denotes the symmetric difference 
of E and B; that is, 

EAB = (E\B) U (B\E) . 

The theorem of Fusco, Maggi and Pratelli [16] says that for all d > 2, there is a constant C(d) 
depending only on d, so that for all Borel sets E in R d , 

f , \ (d-i)/d 

P(E) > [1 + C(d)A 2 (E)]a d IA\E\) . (6.3) 

This improved on an earlier result of Hall, Hayman and Weitsman [191 [20] of the same character, 
but with A 4 (E) in place of A 2 (E). The exponent 2 is sharp; see [26] for further discussion and 
background. 

To apply this result to the shape problem for the local free energy functional, suppose that for 
some trial function m and some e, rj > 0, one has that the volume of C lying outside every ball of 
radius (1 + e)R is at least r)\C\. 

Now let Efi be defined by 

Eh '■= { x : m(x) > h } , 

so that Eh + = C, and Eh_ = A U C. By containment, if k > h, then Ej~ C Eh- This has the 
consequence that for all h G [h~, h + ], the volume of Eh lying outside every ball of radius (1 + e)R 
is at least rj\C\. 

Moreover, since in the critical scaling regime, Lemma 14.61 savs that \A\ is negligible compared 
to \C\, for all sufficiently large L, and all h E [h-,h+], the radius of a ball in M. d with the same 
volume as Eh is no greater than (1 + e)R. 
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Consequently, for each h E A(E; l ) > r\. Therefore, by (|6.3p . and one more use of 

|£ fc | > \E h+ \ = \C\, 

P{E h ) > (1 + C(d)r, 2 )a d R( d - l V d . 

In the context of (II, lip , for almost every h, P{Eh) = \Th\, the d— 1 dimensional Haussdorf measure 
of the set = { x : m(x) = h }, and hence we have a lower bound on \Th\ that is uniform in 
h. Using this in (jl.lip . one gets a lower bound on the free energy of m that is larger than the 
minimizing value by a factor of (essentially) (l + C(d)r] 2 ). Hence if the free energy of a trial function 
is sufficiently close to the minimizing value, r\ must be correspondingly small. This forces C to be 
very nearly a ball. A similar argument applies to each of the for h € (h—,h+), although the 
degree of control on the roundness of E^ diminishes as h approaches h~. A precise statement in 
terms of the LP distance between m and an ideal round droplet profile may be found in |llj . 

It is an interesting open problem to develop bounds of this type that would apply to the Gates- 
Lebowitz-Penrose free energy functional. It would be very surprising if the physical model behind 
it did not capture enough physical reality to control the shape of droplets in near-minimizers. We 
note that one can apply stability for the isoperimetric inequality to the GLP functional, but only 
in the sharp interface scaling limit, as discussed in |30| . However the size of our critical droplet 
goes to zero in this limit. 

What would seem to be useful here would be a stability result for the Riesz rearrangement 
inequality, or at least the special case in which one of the three functions is already rearranged, 
and has all centered balls of sufficiently small radius as level sets. A general stability result for it 
might be quite subtle; the cases of equality were only determined relatively recently by Burchard 
[8]. However, for J is radially symmetric and well-behaved, there is a simpler and cleaner result on 
the cases of equality due to Lieb [22J, and in this setting, one might expect a stability result that 
would force the other two functions to be nearly rearranged for near minimizers. This will be the 
subject of future research. 
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